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INTERMEDIATE AGE STARS AS ORIGIN OF THE LOW VELOCITY DISPERSION NUCLEAR RING IN MRK 1066 
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ABSTRACT 

We report the first two-dimensional stellar population synthesis in the near-infrared of the nuclear region 
of an active galaxy, namely Mrk 1066. We have used integral field spectroscopy with adaptative optics at the 
Gemini North Telescope to map the to map the age distribution of the stellar population in the inner 300 pc at 
a spatial resolution of 35 pc. An old stellar population component (age >5 Gyr) is dominant within the inner 
ss 160pc, which we attribute to the galaxy bulge. Beyond this region, up to the borders of the observation 
field (~300pc), intermediate age components (0.3-0.7Gyr) dominate. We find a spatial correlation between 
this intermediate age component and a partial ring of low stellar velocity dispersions (cr t ). Low-cr t nuclear 
rings have been observed in other active galaxies and our result for Mrk 1066 suggests that they are formed 
by intermediate age stars. This age is consistent with an origin for the low-<x t rings in a past event which 
triggered an inflow of gas and formed stars which still keep the colder kinematics (as compared to that of the 
bulge) of the gas from which they have formed. At the nucleus proper we detect, in addition, two unresolved 
components: a compact infrared source, consistent with an origin in hot dust with mass « 1.9 x 1O _2 M , and 
a blue featureless power-law continuum, which contributes with only *15% of the flux at 2.12/mi. 
Subject headings: galaxies: individual (Mrk 1066) — galaxies: nuclei — galaxies: Seyfert — galaxies: stellar 
content — infrared: galaxies 



1. INTRODUCTION 

Optical spectroscopy on scales of hundreds of parsecs 
around the nucleus of Seyfert galaxies have shown that 
in «40% of them the active g alactic nucleus (AGN) 
and yo ung s tars co-exist (e.g. [Storchi-Bergmann et alj 
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ICid Fernandes et alj 



iGonzalez Delgado. Heckman & Leitheren 
12004b lAsarietal.1 12001 
the so-called 



to 



2001 

Dors et alJ 120081) . providing s upport 
AGN-Starburst connection ( e.g. | Norman & Scovill 
Terlevich. Diaz & Terlevich 1990; Heckman et 



Heckmanj 120041; IRiffel etal] l2009al) . 
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In particular, these 
authors have pointed out that the main difference between 
the stellar population (hereafter SP) of active and non-active 
galaxies is an excess of mainly intermediate age stars in the 
former. 

A similar result has been found in recent SP studies in 
the near-infrared (here after near-IR), using the technique 
of spectral synthesis dRiffel et all l2009bl 12007b . Using 
integrated spectra of the central few hundreds of parsecs 
in a sample of 24 Seyfert galaxies, these authors have 
shown that the continuum is dominated by the contribu- 
tion of intermediate-age stellar population components 
(SPCs). In addition, they found that the near-IR nu- 
clear spectra of about 50 % of the Seyfert 1 and ~20 % 
of the Seyf ert 2 galaxies show emission from hot dust 



or the Seytert / g alaxies show emission from hot dust 
(IRiffel et all l2009bt [Riffel. Storchi-Bergmann & McGregoi 
Riffel et all l2009at iRodrfguez Ardila & Mazzalay 
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2006E IRodrfguez-Ardila. Contini & Viegasll2005afr 

Using a somewhat distinct technique. iDavies et al.l ( 120071) 
obtained near-IR integral field spectroscopy to investigate the 
circumnuclear star formation in 9 nearby Seyfert galaxies at 
spatial resolutions of tens of parsecs. They have modelled 
the Bry equivalent width, supernova rate and mass-to-light 
ratio to quantify the star formation history in the center of 

Electronic address: rogemar@smail.ufsm.br 



these galaxies using their code stars. They found that the 
ages of the stars which contribute most to the near-IR contin- 
uum lie in the range 10-30 Myr, but point out that these ages 
should be considered only as "characteristic", as they have 
not performed a proper spectral synthesis, arguing that there 
may be simultaneously tw o or more SPs that are not coeval 
(IDavies et al.ll27Ml2006l) . 

In this paper we present, for the first time, two-dimensional 
(hereafter 2D) SP synthesis in the near-IR for the inner hun- 
dreds of pc of an active galaxy - namely Mrk 1066 - using in- 
tegral field spectroscopy with adaptive optics, which allowed 
us to derive the contribution of distinct SPCs to the near-IR 
spectra and map their spatial distributions. This paper is or- 
ganized as follows. In Sec. [2] we describe the observations, 
data reduction procedures and the spectral synthesis method; 
in Sec. [3]we present our results, which are discussed in Sec. [4] 
The conclusions are presented in Sec. [5] 

2. THE DATA AND SYNTHESIS CODE 
2.1. Mrk 1066 

Mrk 1066 is an SB0 galaxy harboring a Seyfert 2 nu- 
cleus and located at a distance of 48.6 Mpc, for which 
1 "corresponds to 235 pc at the galaxy. Previous stud- 
ies of the central region of Mrk 1066 have shown that 
line emission from high-excitation gas is dominated by 
gas in a bi-conical outflow oriented along position angle 
PA=135/315°, which seems to be associated to the radio 
jet, while the low-ex citation gas is more restricted to the 



1999; 



2010; 



plane of the galaxy dBower et alj 119951; iNagar et al 
Knop et al.l 1200 ll: IRiffel. Storchi-Bergmann & Nagari 
Riffel & Storchi-Bergmann! 1201 Oh . We selected this galaxy 
for the present study because it shows a signature of the pres- 
ence of young/intermediate age stellar population in its stellar 
kinematics in the inner few hundred parsecs as illustrated in 
the stellar velocity dispersion (cr„) map shown in the top-left 
panel of Fig. [2] This map was obtained from the fit of the 
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CO band heads around 2.3 ami by a combination of stellar 
spectr al templates from IWinge. Riffel & Storchi-Bergmannl 
(2009) using the penalized Pi xel-Fitting (pPXF) method of 
Cappellari & Emsellem (2004). It shows a partial ring of low 
0%, a signature of the presence of stars which still have the 
kinematics of the cold gas from which they have formed 1 . 
Our first goal with the spectral synthesis is to derive the 
age of this population. Additionally, high spatial resolution 
spectroscopy has shown that the near-IR nuclear continuum 
is co nsistent with a large contribution from h ot dust emis- 
sion dRiffel. Storchi-Berg mann & Nagar]|2010|), not d etected 
in lower resolution observations dRiffel et alj 1200 9b). Our 
second goal is to quantify this contribution with spectral syn- 
thesis. 

2.2. Observations and data reduction 

Mrkl066 was observed with the Gemini Near-Infrared 
Integral-Field Spectrograph (NIFS; iMcGregor eTal] 120031) 
operating with Gemini North Adaptive Optics system AL- 
TAIR in September 2008 under the programme GN-2008B- 
Q-30. Two set of observations were obtained, one at the J 
band, covering a spectral region from 1.15 //m to 1 .36 jum with 
a two-pixel spectral resolving power of 6040; and another at 
the Ki band at resolving power of 5290 and covering the spec- 
tral range 2. 1-2.5 fim. 

The data reduction, described in 

iRiffel. Storchi-Bergmann & Nagari (120101) . was performed 
using the gemini iraf package and followed standard pro- 
cedures. The resulting calibrated data cube contains 784 
spectra covering a angular region of 2'.'8x2'.'8 at a sam- 
pling of 07 1x071, which corresponds to 650x650 pc 2 and 
23.5x23.5 pc 2 , respectively. The angular resolution is ^0715, 
corresponding to 35 pc at the galaxy. 

2.3. Analysis 

In order to model the continuum and stellar absorp- 
tion features we used the s tarlight spectral synthesis code 
dCid Fernandes et alJl200l Hool 120091: iMateus et alJ 120061: 
lAsari et al] 120071) . According to the authors starlight mixes 
computational techniques originally developed for semi em- 
pirical population synthesis with ingredients from evolution- 
ary synthesis models. 

In summary, the code fits an observed spectum Ox with a 



ulations ( S SPs) Following IRiffel et al.l d2009bl) we use the 
iMarastonl (12005b Evolutionary Population Synthesis (EPS) 
models as base set (described below). This choice is due to 
the fact that Maraston's models include the effects of stars 
in the TP-AGB phase, wich is crucial to model NIR SP (see 
Riffel Rodrfguez-Ardila & Pastorizdl2006t IRiffel et al.ll2007l 
2008bll2009bh . Extinction is modeled by starlight as due to 
fore ground dust, and parameterized b y the V-band extinction 
A v fcid Fernandes et al.ll2004l 120051) . Note that all the com- 
ponents are reddened by the same a mmount. We use the CCM 
ilCardelli. Clayton & Mathisll 19891) extinction law. 

Basically, the co de solves the following e quation for a 
model spectrum M A dCid Fern andes et al. 20051): 



Mi =M, 



.10 



^Xjb M r A <g> G(v„, cr t ) (1) 

1 A detailed discussion of the stellar kinematics is presented in 
IRiffel & Storchi-Bergmamjl f20TOh . 



were bj^ r A is the reddened spectrum of the y'th SSP normal- 
ized at A ; r A = lcr - 4 ^- 1 ^ 10 ) is the reddening term; M A0 is 
the synthetic flux at the normalisation wavelength; x is the 
population vector; ® denotes the convolution operator and 
G(v„, <x„) is the gaussian distribution used to model the line- 
of-sight stellar motions, it is centred at velocity v s with dis- 
persion o~ *. 

The final fit is carried out with a simulated annealing plus 
Metropol is scheme, which searches for the minimum of the 
equation dCid Fernandes et aTll2005l) : 



X 2 = Yj [(0 ' 1 ~ Mx)wx]2 



(2) 



where emission lines and spurious features are masked out by 
fixing w,i=0. 

The spectral base was constructed wi th EPS models of 
IMarastonl {2005) as in IRiffel et al.l d2009bl) . and comprise sin- 
gle stellar populations (SSPs) synthetic spectra covering 12 
ages (f =0.01,0.03, 0.05, 0.1, 0.3, 0.5, 0.7, 1, 2, 5, 9 and 
13 Gyr) and four metallicities (Z =0.02, 0.5, 1, 2 Z Q ). We 
also included blackbody functions for temperatures in the 
range 700-1400K in steps of 100 K dRiffel et al.ll2009bl) and a 
power-law (F y oc y~ L5 ) in order to account for possible con- 
tributions from dust emission and from a featureless contin- 
uum (FC), respectively, in regions close to the nucleus (e.g. 
ICid Fernandes et alj|2004l) . 

3. RESULTS 

In Fig.Q]we present the results of the synthesis for two spec- 
tra extracted within an aperture of 07 1 x07 1 , one at the nucleus 
(top panel) and the other at 1" south-east of it (bottom panel). 
The observed (black) and synthetic (red) spectra were normal- 
ized to unit at 2.12/im. This wavelength was chosen because 
the K-band spectra present a higher signal-to-noise ratio than 
those in the J-band and the sp ectral region near 2 .12 /mi is free 
of emission/absorption lines (IRiffel et al. 2008b). In the inner 
072 (45 pc) radius the contribution from a FC and from dust 
emission are important, as illustrated in the top panel of Fig.Q] 
where the combined contribution of all SPCs is shown in blue, 
the FC as a dotted line and the composite blackbody function 
(sum over all temperatures) as a dashed line. We point out that 
the syntetic spectra (sum of all components) was reddened by 
the reddenin g amount detected in the ob served spectra. 

Following lCid Fernandes et alj d2004l) . we have binned the 
contribution of the SPCs (xj) to the flux at 2.2yum into four 
age (f) ranges: young (x v : t < 100 Myr); young-intermediate 
(x yi : 0.3 < t < 0.7 Gyr), "intermediate-old (x io : 1 < t < 2 Gyr) 
and old (x a : 5 < t < 13 Gyr). The percent flux contribution at 
2.12jum from stars of each age range are shown in Fig. [2] on 
average, the young-intermediate SPCs dominate the contin- 
uum emission from the central region of Mrk 1066, but there 
is significant variation over the NIFS field. While in regions 
farther than r ~077 - 160pc from the nucleus - the contri- 
bution of the young-intermediate age stars reaches values of 
up to 100%, closer to the nucleus its contribution is negligi- 
ble. Within this region the stellar population is dominated by 
the old component, whose contribution reaches « 50 %, fol- 
lowed by the young component, whose contribution reaches 
~30 % of the flux at 2. 12yum. Also within this region, but un- 
resolved by our observations (thus within 0715 from the nu- 
cleus) we find a contribution of up to 50% from the combined 
blackbody components. Finally, we have found a small con- 
tribution from the AGN FC component, which reaches at most 
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X(jj,m) 

Fig. 1. — Sample fits of the SPCs to the nuclear spectrum (top) and to an extra-nuclear one at l'.'O south-east of the nucleus (bottom). The observed spectra are 
shown in black and the fits in red. In the top panel the dashed line shows the contribution of hot dust, the dotted line that of the featureless continuum and the 
blue continuous line shows the contribution of the combined SPCs. 



15 % at the nucleus, also unresolved by our observations. The 
intermediate-old component contributes with less than 20 % 
at most locations, with larger values more concentrated to the 
west-northwest of the nucleus. 

Besides the SPC distributions, the synthesis outputs the 
average reddening of the stellar populations, shown in the 
central panel of Fig. [3] It reaches the highest values, of 
up to E(B - V) = 1.4, along the position angle 128°, 
which is the orientation of the li ne of nodes of the galaxy 
dRiffel & Storchi-Bergmannl2010l) . 

The robustness of the SP fit can be measured by the 
percent mean deviation (adev), \0,\ - M,\\/0,\, where Ox 
is the observed spectrum and Mx is the fitted model 
(ICid Fernandes et all 12004 120051) . The resulting adev map 
is shown in the right panel of Fig. [3] and presents values 
adev < 2.5 % at most locations, indicating that the model 
reproduces very well the observed spectra. Nevertheless, at 
regions close to the border of the NIFS field the adev reaches 
values of up to 5%. At these locations we found an increase 
in the contribution of the young component x y , which has 
thus larger uncertainties than in the rest of the NIFS field and 
should be further investigated with better data. 

4. DISCUSSION 



In general, our results are in good agreement 
with those of previous near-IR studies fo r sin- 
gle aperture nuclear spectra dRiffel et al.1 l2009bt 
iRamos Almeida. Perez Garcia & Acosta-Pulidol l2009t) . 
The main novelty of our work is the 2D mapping of the 
stellar population using near-IR spectra, which shows spatial 
variations in the contribution of the SPCs in the inner few 
hundred parsecs of Mrk 1066. 

The significance of the above variations is further enhanced 
by the comparison with the <x„ map presented in the top-left 
panel of Fig. [2] This map shows a partial ring of low cr* 
values (» 50kms~') surrounding the nuclear region (which 
has <x„ a 100 km s _1 ) at « 1" from it. Such rings are 



com monly observed in t he central region of active galax- 
ies dBarbosa et all ]2006t iDeo. Crenshaw & Kraemerl 120061 : 
ISimoes Lopes et alj|2007t iRiffel et al.ll2008al l2009al) and in- 
terpreted as being due to colder regions with more recent 
star formation than the underlying bulge. The comparison of 
the results from the SP synthesis of Mrk 1066 with the cr, 
map shows that the low cr,-values ring is associated with the 
young-intermediate age SPC, while the highest cr„ values are 
associated with the old SPC, supporting the use of low stellar 
velocity dispersion as a tracer of younger stars in the bulge 
and confirming the interpretation of the above studies. 

The flux-weighted SPC contributions depend on the nor- 
malization point and thus the comparison with results from 
other spectral regions should be done with caution. A phys- 
ical parameter which does not depend on the normaliza- 
tion point and spectral range used in the synthesis is the 
mass of the stars. Thus, we constructed maps for the mass- 
weighted contribution of each SPC. The contributions for the 
young population (m y ), young-intermediate population (/«,,), 
intermediate-old population (m J0 ) and old population (m ), 
are shown in the rightmost panel of second row, bottom-left, 
bottom-central and bottom-right panels of Fig. [2] respectively. 
The mass-weighted contribution of the young population is 
very small over the whole field, while m a dominates within 
~160pc from the nucleus and m v , dominates in the circumnu- 
clear ring. 

The only previous 2D stellar population studies of active 
galaxies in the near-IR published to date are those fr om the 
group of R. I. Davies. In particular. lDavies et al.l d2007) inves- 
tigated the circumnuclear star formation in nine Seyfert galax- 
ies using near-IR IFU observations and also found circumnu- 
clear disks of typical diameters of tens of pc with a "char- 
acteristic" age in the range 10-300 Myr based on measure- 
mentes of the Bry emission-line equivalent width, supernova 
rate and mass-to-light ratio. Thus, the results we have found 
for Mrk 1066 are in reasonable agreement with those found by 
Davies group. Nevertheless, the methodology adopted in the 
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Fig. 2. — Stellar velocity dispersion map (top-left panel) and spatial distributions of the percent contribution of each SPC to the flux at A = 2.12/jm (xj) and to 
the mass (nij), where j represents the age of the SPC: young (y: < 100 Myr), young-intermediate (yi: 0.3-0.7 Gyr), intermediate-old (io: 1-2 Gyr) and old (o: 
5-13 Gyr). 
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Fig. 3. — From left to right: Spatial distribution of the dust emission summed over all temperatures, reddening map and adev map (percent mean deviation from 
the spectral fit). 
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present work allowed us not only to obtain a "characteristic" 
age, but to also map, for the first time, the spatial distribu- 
tion of stars of different ages in the central region of a Seyfert 
galaxy, on the basis of n ear-I R integral-field spectroscopy. 

As presented in Sec. 12.31 our spectral base include four 
distinc metalicities: Z_=0.02, 0.5, 1, 2 Z . Following 
ICid Fernandes et ail (120041) . we can estimate the mean flux 
and mass-weighted metalicity in the fit by 

N 

</. >,-X^ ( 3 ) 

7=1 

and 

if, 

< Z„ > M = J] m ! Z h ( 4 ) 
7=1 

where < Z„ >i is the flux-weighted mean metalicity, < Z* >m 
is the mass-weighted mean metalicity and N* is the num- 
ber of SSP in the spectral base. In the case of Mrk 1066 
we obtain a nearly solar metalicity - < Z* >i- 0.021 and 
< Z„ >m— 0.017 - in good agreement with those found by 
iRiffel et all (12009b) from near-IR spectral synthesis using a 
single aperture nuclear spectrum. 

In further support of the results of the synthesis, we found 
that the average reddening map derived for the stellar popula- 
tion (central panel of Fig. [3} is in close agreement to the one 
derived for the narrow-line region using emission-line ratios 
( IRiffel. Storchi-Bergmann & Nagarll2010l) . presenting a simi- 
lar S-shaped structure. 

The synthesis also confirmed the presence of an unresolved 
blackbo dy (BB) component at the nucleus , whic h we had 
found in IRiffel. Storchi-Bergmann & Nagari (120101) from the 
fit of the nuclear spectrum by a BB function plus a power-law. 
On the ba sis of its luminosity , we have followed the calcula- 
tions as in Riffel et al. (2009b) to estimate a total dust mass of 
M HD ~ 1-9 x 10' 2 M R , which is in good agreement with the 
value obtained in IRiffel. Storchi-Bergmann & Nagari (120101) 
and with t h ose ob se rved for other Seyfert galaxies (e.g. 
Riffel et all l2009at IRiffel. Storchi-Bergmann & McGregor 
20091: IRiffel etalj l2009bt iRodrfguez Ardila & Mazzalav 
20061 iRodrfguez-Ardila. Contini & Viegasll2005al) . 

As discussed in Sec. [3] we found that the FC com- 
ponent contributes with ~ 15% of the 2.12yum nuclear 
continuu m, which is in goo d agreement with the result 
found by Riff el et ail J2009b). This contribution is never- 
theless uncertain, as it is hard to distinguish a reddened 
young st arburst (t < 5Myr) from an AGN FC compo- 
nent (e.g. IStorchi-Bergmann et al.ll2000t ICid Fernandes et all 
I2004t IRiffel et al.l l2009bl) . On the other hand, the con- 
tribution of the young population x y is dominated by 
stars with ages of 50-100Myr, thus there does not seem 
to be a significant contribution from very young stel- 
lar populations, suggesting that the FC component orig- 
inates from the emission of the AGN. This interpreta- 
tion is supported by the detection of weak broad com- 
ponents for the Pa/3 an d Bry emission lines for the nu- 
clear s pectra of Mrk 1066 ([Riffel. Storchi-Bergmann & Nagari 
12010b IVeilleux.Goodrich & Hilll 119971) . However, spectro- 
polarimetric studies, show n o evidence for the hidden Seyfert 
1 nucleus in Mrk 1066 (e.g. lBian & Gull2007l) . Thus, the ori- 
gin of the FC in Mrk 1066 must be further investigated. 

We finally point out one caveat of the stellar population syn- 
thesis in the near-IR. In IRiffel. Storchi-Bergmann & Nagaii 



(120101) . on the basis of emission-line ratios, we have iden- 
tified a large star-forming region at (y.'5 south-east of the 
nucleus. In order to be active, this region should have an 
age Z IP Mvr (IRiffel et al.ll2009atlDo7s et al.l l200a[ Dfaz et alj 
l2007t iKennicutt et al.ll 1989^ 7 Nevertheless, our synthesis did 
not find a young stellar component at this location, indicat- 
ing that spectral synthesis in the near-IR (at least in the J and 
K bands) is not a good tracer of very recent star formation 
(age < lOMyr). This result is expected since the young (blue) 
population has its peak emission in ultraviolet/optical wave- 
lengths, while in the near-IR its contribution to the flux is very 
small, being hard to detect this component in this wavelength 
region. 

In order to have a complete census of the stellar population 
components it is thus essential to study the stellar content of 
the central region of active galaxies using distinct methods 
and spectral regions. 

5. CONCLUSIONS 

The present work reports for the first time spectral synthe- 
sis in the near-IR with 2D coverage for the nuclear region of a 
Seyfert galaxy (Mrk 1066) within the inner a 300 pc at a spa- 
tial resolution of «35pc. We have mapped the distribution 
of stellar population components of different ages and of their 
average reddening. The main conclusions of this work are: 

• The age of the dominant stellar population presents 
spatial variations: the flux and mass contributions 
within the inner « 160 pc are dominated by old stars 
(f >5 Gyr), while intermediate age stars (0.3 < t < 
0.7 Gyr) dominate in the circumnuclear region; 

• There is a spatial correlation between the distribution of 
the intermediate age component and low stellar velocity 
dispersion values which delineate a partial ring around 
the nucleus. Similar structures have been found around 
other active nuclei, and our result for Mrk 1066 sug- 
gests that these nuclear rings (and in some cases disks) 
are formed by intermediate age stars. 

• There is an unresolved dusty structure at the nucleus 
with mass Mhd * 1.9 X 10~ 2 M , which may be the 
hottest part of the dusty torus postulated by the unified 
model of AGN and a small contribution from a power- 
law continuum («T5 % of the flux at 2.12//m); 

• The near-IR synthesis seems not to be sensitive to very 
recent star formation (with t < 5 Myr), reinforcing the 
importance of multi-wavelength stellar population stud- 
ies of the central region of active galaxies. 



Based on observations obtained at the Gemini Observa- 
tory, which is operated by the Association of Universi- 
ties for Research in Astronomy, Inc., under a cooperative 
agreement with the NSF on behalf of the Gemini partner- 
ship: the National Science Foundation (United States), the 
Science and Technology Facilities Council (United King- 
dom), the National Research Council (Canada), CONICYT 
(Chile), the Australian Research Council (Australia), Min- 
isterio da Ciencia e Tecnologia (Brazil) and south-eastCYT 
(Argentina). This work has been partially supported by the 
Brazilian institutions CNPq and CAPES. 
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